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The so-called "Quantum Inequalities" , and the "Quantum Interest Conjecture" , use quantum field 
theory to impose significant restrictions on the temporal distribution of the energy density measured 
by a time-like observer, potentially preventing the existence of exotic phenomena such as "Alcubierre 
warp-drives" or "traversable wormholes" . Both the quantum inequalities and the quantum interest 
conjecture can be reduced to statements concerning the existence or non-existence of bound states 
for a certain one-dimensional quantum mechanical pseudo-Hamiltonian. Using this approach, we 
shall provide a simple variational proof of one version of the Quantum Interest Conjecture in (3+1) 
dimensional Minkowski space. 
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I. INTRODUCTION 

It is "well-known that quantum physics permits arbi- 
trarily large negative energy densities at individual points 
0, 0] , though averages and total energies over volumes or 
lines are much more tightly constrained. This is of crit- 
ical importance when developing singularity theorems, 
and other theorems based on global analysis, in that this 
simple observation is enough to guarantee that the so- 
called "classical energy conditions" are not fundamental 
physics; they are at best classical approximations to a 
more subtle quantum universe Q. 

Note that without something similar to the energy con- 
ditions to constrain the spacetimes one wishes to consider 
as "physical" , one can construct arbitrarily weird space- 
times containing such exotic objects as warp-drives [J, Q , 
traversable wormholes Qj H, H, 0| , singularity- free "black 
holes" [l0[ , GNACHOs |nf. violations of the generalized 
second law [l2[ , violations of cosmic censorship [l3| , and 
even time machines [3, [lfl [l6j]. Because of the need 
to keep such oddities somewhat constrained, and avoid a 
complete free-for-all, much work has gone into develo pin g 
"averaged energy conditions" of various types [l^lillll9j . 
though even here there are issues associated with quan- 
tum field theory anomalies [2(| ■ 

In particular, within the framework of semi-classical 
General Relativity (GR), the Quantum Inequalities (QIs) 
are an extremely important tool for constraining both 
negative energies, and more generally exotic phenomena 
which violate the classical energy conditions. Originally 
introduced by Ford [21], and subsequently extensively in- 
vestigated by Ford and Roman [22|,[23|], and their collabo- 
rators 0, [HI , the QIs impose a lower bound for the inte- 
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gral of the renormalized stress-energy tensor along a suit- 
able time interval, weighted by a suitable test function. 
Initially the test function was chosen to be a Lorentzian 
function [23| . but subsequently Flanagan found a result 
independent of the specific choice of the test function 
in flat (1+1) space-time [26j |, With a distinct approach, 
Fewster and Eveson also obtained a bound independent 
of the weight function in 4 and 2 dimensional Minkowski 
space [2^ , though their result in 2 dimensions was some- 
what weaker than that of Flanagan [26| . 

A related but distinct restriction for negative energies 
is imposed by the Quantum Interest Conjecture (QIC), 
introduced by Ford and Roman [28]. Informally speak- 
ing, the QIC states that "overall" the energy density 
must be positive, any negative energy density in one 
region must be overcompensated by positive energy in 
another region. When restricted to isolated pulses of en- 
ergy, the QIC not only constrains the amount of negative 
energy in a pulse, but also the time interval between any 
negative energy pulse and the larger positive energy pulse 
that must also be present. The (positive) nett energy of 
the two pulses is called the Quantum Interest, and grows 
monotonically with the time separation. Initially Ford 
and Roman used delta function energy pulses to prove the 
conjecture, and it had not been generalized to arbitrar- 
ily shaped energy pulses in 4 dimensional flat space-time, 
nor was it initially formulated for curved space-time. 

A much more general approach to the QIC that is not 
restricted to isolated pulses was developed by Fewster 
and Teo in [2{J, wherein they created a Sobolev-space- 
based technical apparatus to reinterpret the QIs and the 
QIC as eigenvalue problems for a one-dimensional quan- 
tum mechanical pseudo-Hamiltonian. This technique re- 
lates the (1+1) dimensional QIs and QIC to the nonex- 
istence bound states for the Schrodinger equation (SDE) 
in one dimension, while the (3+1) dimensional QIs and 
QIC are related to the nonexistence bound states for 
the bi-harmonic Schrodinger equation (bSDE). Within 
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this framework, Fewster and Teo [29| , and later Teo and 
Wong [30f , were able to prove the QIC in (1+1) dimen- 
sional Minkowski space; however they could not extend 
their proof to the (3+1) dimensional case. In Teo and 
Wong's approach [301 ] critical use is made of a theorem 
by Simon [31j , regarding the existence of bound states for 
the 1 dimensional SDE. That is, Simon's theorem is used 
to prove a version of the QIC in (1+1) dimensions, but 
the generalization to (3+1) dimensions is troublesome. 

Furthermore, in [29j Fewster and Teo proved that is not 
possible for a positive delta function pulse to compensate 
a negative delta function pulse in 4 dimensions, regard- 
less how small the time separation between them could 
be. Additionally, in [32| Pretorius found an upper bound 
for the separation between two general energy pulses by 
applying a scaling argument to the test function. With 
this strategy, Pretorius also showed that the Quantum 
Interest grows almost linearly as the pulse separation in- 
creases, which is the same result obtained in [281 ] . 

Turning to the current article: Using the approach ex- 
plained in [2^ | , we shall provide a simple variational proof 
for a variant of the QIC in 4 dimensional flat space- 
time, (and in fact in any even dimensional spacetime), 
by proving the equivalent of Simon's theorem for the bi- 
harmonic Schrodinger equation (and more generally for 
the multi-harmonic Schrodinger equation). We shall do 
this via a variational argument coupled with a power se- 
ries expansion of the test function used in the related 
one-dimensional pseudo-Hamiltonian problem. 

The layout is as follows: In section II we explore 
the QIs and the QIC from the point of view of eigen- 
value problems, to obtain the related Schrodinger equa- 
tions (SDE and bSDE). Then in section III we will use 
an appropriate class of test functions to prove Simon's 
theorem via variational techniques applied to a power 
series expansion, and then in section IV repeat this 
strategy for the (3+1) dimensional related problem to 
prove the equivalent of Simon's theorem for the bSDE. 
(Furthermore, as shown in the appendix, this naturally 
leads to a version of Simon's theorem for multi-harmonic 
Schrodinger equations.) When translated back to the 
underlying even-dimensional Minkowski spacetime this 
guarantees the existence of quantum interest, though this 
particular calculation is mute as to the amount of interest 
and "interest rate" . We conclude with a brief discussion 
in section V. 



II. QUANTUM INEQUALITIES AND THE 
QUANTUM INTEREST CONJECTURE 

A. The Quantum Inequalities 

The Quantum Inequalities give a lower bound for the 
expectation value of the renormalized stress-energy ten- 
sor (in a quantum state %p) when evaluated along a com- 



plete timelike geodesic, 

poo 

h,« = I (T^ n (t,0))^w(t)dt (1) 

Xj 

\T^{t^))^w{t)dt. (2) 

Here w(t) is a more-or-less freely-specifiable test func- 
tion (weighting function), which is non-negative and in- 
tegrates to unity: w{x) dx = 1. Initially, the stress- 
tensor is treated in the test- field limit, so the background 
is simply taken to be Minkowski space, and the timelike 
geodesic is specified by the 4- velocity V a — (1;0). Since 
our integrals almost always run from —00 to +00 we shall 
for brevity often use the shorthand symbol § = J^°° ■ 

In their initial analysis, Ford and Roman used a 
Lorentzian test function to yield to a lower bound which 
depend on modified Bessel functions [23j. However a 
more general inequality was found by Flanagan in (1+1) 
dimensional flat space-time [26| , where he used Quantum 
Field Theory for a massless scalar field minimally coupled 
to gravity to find the optimum lower bound independent 
of the specific choice of the test function, 
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K(t)) 2 

w(t) 



dt. 



(3) 



Fewster and Eveson used a distinct approach [27[ to yield 
inequalities in (1+1) dimensions, 



1 r {w'(t)) 2 , 

and in (3+1) dimensions, 



■ib/(^w) ! 



dt. 



(4) 



(5) 



It is obvious that the better bound in (1+1) dimensions 
is given by Flanagan's result rather than the one by Few- 
ster and Eveson. However the most general and optimum 
bound in (3+1) dimensional flat space-time, for a mass- 
less scalar field minimally coupled to gravity, is shown in 
equation ([5]). 

Combining the results of Flanagan, and Fewster and 
Eveson, is possible to rewrite the QIs in 2m-dimensional 
spacetime in the more general form [29| . (writing w(t) — 
fit) 2 to automatically enforce the positivity constraint), 



W(t,o)L [f(t)] 2 dt>-— i[D m f{t)fdt, (6) 



where D is the derivative operator, and the various dif- 
ferent constants are given by 
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m7r m-l/2 2 2m r ( m _ 1 



m = 1; 
) m > 2; 



(7) 



and where the weighting function is now normalized to 
J_°° \f{x)\ 2 dx = 1. The case m = 1 is equivalent to 
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equation ([3]), while equation ((5|) is equivalent to m = 2. 
After integrating equation © by parts, the QIs become a 
statement regarding the non-existence of negative eigen- 
values for a one-dimensional pseudo-Hamiltonian, 



where, 



W\f) > 0, (8) 



H=(-irD 2m +Cm (nr (t,o)) v ,. (9) 



The quantity c m ( T[^ n (t, 0) )^ can effectively be viewed 
as an "potential" V for a quantum mechanical system, 

v^c m (r —(t,o))^, (io) 

and the QIs become the statement 

{f\{(-l) m D 2m + V}\f)>0, (11) 

We shall now use the framework of quantum mechanics 
to analyze the QIs and the QIC as an eigenvalue problem. 



B. The QIC as an eigenvalue problem 

As adopting this point of view allows us to use the 
mathematical background of one-dimensional quantum 
mechanics, we will change to the standard quantum me- 
chanical notation. For instance, the differential operator 
([9|) becomes the quantum pseudo-Hamiltonian 



H = P 



V, 



(12) 



where P and V are operators on the usual Hilbert space 
of square- integrable functions. The eigenvalue prob- 
lem for this Hamiltonian, in coordinates, is the ODE 
[ordinary differential equation] (we have effectively set 
h/2m — > 1 to simplify the algebra), 



(-iy 



dx 2r ' 



ip(x) + V(x) <p(x) = Eip(x), 



(13) 



for the eigenfunctions ip(x) which again belong to the 
usual Hilbert space of square-integrable functions. (If 
desired, a more rigorous approach can be followed, as 
in HH). Taking m — 1, (i.e., (1+1) spacetime dimen- 
sions) , we recover the one-dimensional time- independent 
Schrodinger equation (SDE) 



dx 2 



<p(x) + V(x) (p(x) = E ip(x), 



(14) 



while, for m — 2, (i.e., (3+1) spacetime dimensions), 
we obtain the one-dimensional time-independent bi- 
harmonic Schrodinger Equation (bSDE), 



-j^j ip(x) + V{x) tp(x) = E <p(x). 



(15) 



For the case of the SDE, there is a theorem by Simon [3l| 
which, by imposing a constraint on the potential, ensures 



the existence of a bound state for the Hamiltonian. In a 
simple form it reads: 

Simon's Theorem: 

Let V(x) obey (1 + x 2 ) \V{x)\ dx < oo, with V{x) 

not zero almost everywhere. Then H = — d 2 /dx 2 + V(x) 
has a negative eigenvalue if 



V(x) dx < 0. 



(16) 



□ 



With the help of this theorem, it is possible to generalize 
the QIC to a wider set of energy pulses, and not only be 
restricted to 5-function pulses as in [28| . 

In order to better understand the QIC in this frame- 
work, first we have to point out that a potential which 
fulfills the QIs must violate condition (fTB|) . Now let us 
rewrite the potential as 



V(x) = V{x)+ - V(x) — 



(17) 



This splits the potential into its positive part minus its 
negative contributions, with V(x)± > 0. Then, if the 
potential violates condition (|16p . from (|17[) we have 



V{x)+ dx > & V{x)_ dx. (18) 



This means, if we go back to the GR terminology, that 
the QIs imply that the expectation value of the renor- 
malized stress-energy tensor must satisfy 



(19) 



which is slightly stronger than the averaged weak energy 
condition (AWEC). Splitting the energy density into its 
positive and negative parts 



/ rpren t 
\ J 00 



l M); 

this implies 



' rpren 
> J 00 



(*,o): 



'T™{t,Q))_, 



(20) 



(TS§ n (t,0)) + dx>6 (T — (t,0))_ dx. (21) 



That is, the positive part of a stress tensor that satisfies 
the QIs, must always overcompensate its negative part, 
this is one version of QIC. (Of course, there are signifi- 
cant elements missing from this form of the conjecture, 
as compared with the original form and Pretorius' scal- 
ing argument [32j . Specifically there is no immediate way 
to extract a bound for the time separation between the 
energy pulses and the Quantum Interest as a function of 
this separation. However, a benefit off the current dis- 
cussion is that we are not limited to delta function pulses 
of pulses of compact support.) 

In the next section we will prove Simon's theorem for 
the SDE by expanding an appropriate class of test func- 
tions in a power series; this strategy will be used again 
for the bSDE in section 4, and in the appendix a simi- 
lar result will be obtained in arbitrary even-dimensional 
Minkowski spacetime. 
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III. 1+1 DIMENSIONAL MINKOWSKI SPACE 



B. Generic class of test functions 



In order to prove Simon's theorem, avoiding most of 
the technical issues arising in the formal proof [3l[, and 
keeping in mind the idea of generalizing the theorem to 
the bSDE, and even higher-derivative "Hamiltonians" , 
we shall use a power series expansion of the test function 
in the SDE to find the conditions for the potential V to 
bind. 



A. Gaussian wave- function 

Let us start, for simplicity and clarity, with a Gaussian 
test function, 



<Ptest — 



exp{-(x- At ) 2 /(2a 2 )} 



27TCT 



1/2 



(22) 



which automatically enforces the unit normalization 
§ \ip{x)\ 2 dx = 1. Therefore, a variational argument im- 
plies that for the lowest eigenvalue of the SDE one has: 



E < d> [ip'(x) 2 + V(x) ip{x) 2 } dx. 



(23) 



Then the kinetic term is <f>(p'(x) 2 = 1/er 2 , whereas (as- 
suming all the moments | J> x 2n V(x) dx\ < oo so that the 
expansion makes sense) 



V(x) (p(x) 2 dx = 

1 (-1)" 



2na ^ n!(2a 2 )" 

n— 



x <p V{x) {x - nY n dx. (24) 



Interchanging the integral and the summation is justified 
since the Gaussian function ^p(x) 2 is an analytic function 
of x whose Taylor series has an infinite radius of conver- 
gence. 

For a fixed fi we now see that: 



Ea l < 1 



V(x) dx + 0(l/a) 



From this last expression, we can see that if 



V(x) dx < 0, 



(26) 



then for a sufficiently large we can guarantee E a 2 < 0, 
(hence E < and the potential will bind, and so vio- 
late the QIs). This is close to Simon's condition. (We 
currently have a <, rather than Simon's <. The only 
"difficult" thing about Simon's theorem is proving the 
existence of a bound state in the marginal = case. Note 
that the weaker result we have here is already enough, 
in GR language, to assure that the stress tensor satisfies 
the AWEC.) 



To deal with the marginal case, we proceed by using 
the fact that we are free to choose the test function in 
a more-or-less arbitrary manner, and so to use this free- 
dom obtain tighter constraints. Let us start by picking 
some function g{x) that is analytic on the entire real line 
(infinite radius of convergence) and then constructing the 
piecewise analytic function 



h(x)=g(\x\). 



(27) 



Note that h(x) need not, and typically will not, be ana- 
lytic at zero. We then enforce the normalization 



h(x) dx = 1, fx h(x) dx — 0, f x 2 h{x) dx — 1. 



Suitable examples of such functions are 
exp(-V2>|) 



h(x) 



V2 



and 



N 



K x ) = 5Z Ci ex p( _ MM)' 



(28) 
(29) 

(30) 



subject to 



N 



N 



(31) 



We now choose 



h([x- fj]/a) g{\x- n\/a) 

ftest = V = V , (32) 



whence 



(25) and 



\<p(x)\ 2 dx = 1, & x \ip{x)\ 2 dx = fi, (33) 



(34) 



1 (x — fi) 2 \ip(x)\ 2 dx = a 2 . 

Moreover, § ip'(x) 2 dx = n/a 2 , with 

1 t h'{x) 2 , 
' dx. 



4 / h(x) 



(35) 



The numeric value of k depends on the choice of the test 
function and further optimizations may be useful. For 
now, let us focus on the Taylor expansion of g(x): 



n=0 



(36) 
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This is assumed to exist (by analyticity) , and to converge 
on the entire real line. We now deduce the existence of a 
similar power series expansion for h(x): 



h{x) = g(\x\ ) = ^o„ 



(37) 



n=0 



This is now an expansion in the variable \x\. Since it 
depends on \x\, this is not a Taylor series for h(x), but 
by construction it is convergent over the entire real line. 
Then for the shifted and rescaled test function we have 



Wtest? = ~ 

a 



■ Ml 



(38) 



Assuming all the moments | § \x\ n V(x) dx\ < oo, so that 
the expansion makes sense, we furthermore have 



oo 

„ K 1 v— ^ a 

I - - • - 

(T 



i=0 



dx. 



(39) 



In interchanging the summation and the limit we have 
first split § = = + f + °° , and then appealed 

to the analyticity of g(x) and g(— x) on the appropriate 
ranges, finally recombining the two sets of integrals to 
run over the entire real line. 

We now have the potential for additional relevant 
terms to appear in the series expansion for the variational 
bound on the lowest eigenvalue E, and is convenient to 
place some restrictions on the coefficients a n appearing 
in the series expansion. 

First, if ao = 0, then 



h(x) ~ ax \x\ + 



so that 



h(x) 



+ 



(40) 



(41) 



Thus in order for the integral defining k [equation (|35[) ] 
to converge at x = it is necessary to chose a\ = 0. So 
in fact 



and 



h(x) ~ 02 x 2 + 

[h'(x)} 2 



h{x) 



4a 2 



(42) 



(43) 



In this situation we then cannot extract much informa- 
tion, since the SDE now merely gives 



Eg 1 < n + Oio-- 1 ). 



(44) 



Hence, to extract useful information, we need ao =/= 0, 
and so (since probability densities are always positive), 
a Q > 0. 



In contrast, for ao > 0, the integral defining k [equation 
([33])] always converges at x = 0, since 



h(x) ~ ao + ai |x| H , 



and so 



Thus we have 



WW? . aj 
h(x) a 



(45) 



(46) 



Ea 1 < n+a 



a o j> V{ x ) dx+ai <j> \x—fi\ V(x) dx+0(l/(x). 

Therefore, with a sufficiently large the only important 
contributions come from the first two terms (n = 0, and 
n = 1) of the power series expansion. 

Here, as for the Gaussian test functions, the condition 
§ V(x) dx < implies binding; which means that the ex- 
pectation value of the renormalized stress-energy tensor 
must fulfill the AWEC. However, if we consider the bor- 
derline of AWEC violation, i.e., § V(x) dx — 0, then the 
next term in the expansion is now no longer neglectable: 

Ea 2 < K + ax I \x - fi\V(x) dx + 0(l/a). (48) 

But recall that completely a\ is arbitrary, both in sign 
and in magnitude. That is, if for any value of jx 



■fi\V(x) dx J= 0, 



(49) 



then this implies the existence of a bound state for the 
SDE. This is a very strong constraint on the potential and 
it quickly yields a proof of Simon's theorem. Taking the 
converse of the above, if § V(x) dx = then to prevent 
the occurrence of a bound state we must have 



V/x: j \x - fi\ V(x) dx = 0. (50) 
Differentiating twice with respect to /i we get 



V/x : 



2 <t> S(x - n) V{x) dx 



(51) 



Thus any non-zero potential binds if its integral is null. 
Combining this with the fact that § V{x) dx < also im- 
plies binding, now provides a simple variational proof of 
Simon's theorem, equation (|16[) . though under the much 
stronger technical conditions that all the moments exist 
| § \x\ n V(x) dx\ < oo. It is these stronger technical con- 
ditions — which are still physically quite reasonable and 
are certainly satisfied by isolated pulses of stress energy 
— that will make it easy for us to extend Simon's the- 
orem to the bi-harmonic Schrodinger equation, and so 
lead to a (3+1) dimensional version of the QIC. 
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IV. 3+1 DIMENSIONAL MINKOWSKI SPACE 



function, as we used for the SDE, (l32l). Note 



To generalize the condition given by Simon's theorem 
to the bSDE, we proceed as before, expanding appropri- 
ate test functions in power series. Applying a variational 
argument to the bSDE we find that the lowest eigenvalue 
satisfies 



E< d> [ip"(x) 2 +V(x) <p(xf]dx, 



(52) 



assuming that the test functions are normalized. We now 
use this relation to probe for the existence of a bound 
state for the pseudo-Hamiltonian for the bSDE, which is 
ultimately related to the 4 dimensional flat space-time 
QIs and the QIC. 



A. Gaussian wave- function 

We start with the normalized Gaussian test function 
previously used for the SDE, equation (|22|) . The kine- 
matic term yields § ip"(x) 2 dx — k 2 /ct 4 , with k a numeric 
constant. Assuming all appropriate moments exist, we 
again expand the test function in a power series, whence 



V(x) (p 2 (x) dx = 

1 00 (—l\ n 

Eii x <t>V(x)(x-») 2n dx. (53) 



n— v ' 



Then, from (|52j) 
Ea 4 < k 2 + 



2V2vr 



V(x) dx 



V(x) (x- pt) 2 dx + 0(l/a). (54) 



As for the SDE, in this case also § V(x) dx < implies 
that the potential binds for a sufficiently large a. This 
gives us most of Simon's theorem, apart from the ex- 
tremal case where the integral vanishes. 

However, even with the Gaussian test function, in this 
bSDE case we recover more information than for the 
SDE. Specifically if § V (x) dx = 0, then the next term of 
the expansion is important, and to guarantee the absence 
of a bound state we must also enforce: 



V/x : 



(x - y) 2 V(x) dx < 0. 



(55) 



This constraint, while certainly significant, is not quite 
strong enough to imply Simon's theorem. 



B. Generic class of test functions 

Finally, for the bSDE, we consider the same normal- 
ization, (|28p. and the same choice for the generic test 



ip"(x) dx = -j, 

(7 



with 



(x) > dx 



1 

To 



Vh 



[2h{x)h"{x) - h'(xf 
h{xf 



■ dx. 



(56) 



(57) 



(58) 



In view of the fact that we have chosen h(x) = g( \x\ ) 
this becomes (note the delta function arising from twice 
differentiating the absolute value) 



1 

= — x 
16 

{2g(\x\) [g"(\x\)+g'(\x\)5(x)]~g'(\x\) 2 Y 
2(M) 3 



(59) 



dx. 



In order to derive a useful bound we will want k to be 
finite, which means that we want the coefficient of the 
delta function to be zero, that is a% — 0. To make the 
integral converge at zero we also want ao ^ 0, and in 
fact must have ao > 0. In contrast, a 2 and 03 are uncon- 
strained as to sign and magnitude and we have 



Eg 4 < k + a Q a 3 <j> V(x) dx 

+ a 2 a f V(x) |x — fi\ 2 dx 



V(x) \x - fi\ 3 dx + 0(1/ a). (60) 



As before, by suitably choosing a sufficiently large a, we 
see that § V(x) dx < implies binding. 

To extract further information from (I60|) it is neces- 
sary to check what happens if we set § V(x) dx = 0. 
Thereafter the next two terms in the expansion become 
relevant, however both the sign and magnitude of ai and 
a-$ are freely specifiable, hence (for any jj) either 



|x - m| 2 V(x) dx ^ 0, 



|x - fx\ 3 V(x) dx ^ 0, 



(61) 



(62) 



is a sufficient condition for the Hamiltonian to have a 
bound state. 

Conversely, if § V(x) dx — 0, a necessary condition for 
the absence of a bound state is that 

V/i : j> \x - n\ 2 V(x) dx = 0, (63) 

and 

Va* : j> \x - fi\ 3 V(x) dx = 0. (64) 
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Now if we differentiate the last expression with respect 
to (i then 

V,: f\x-rf S i g n(x-»)V(x)dx = 0) (65) 
whence, combining with Eq. (f63|) , we have 

poo 

V^: / \x - fi\ 2 V(x) dx = 0. (66) 
Jfj, 

If we repeatedly differentiate the last expression with re- 
spect to \x then 

V^: / \x - fi\ V(x) dx = 0. (67) 

J (J. 

V/z : / V(a;) da; = 0. (68) 
J ii 

and finally V(x) = 0. That is, if /V(ac) da; = 0, a 
necessary condition for the absence of a bound state is 
that V(i) = 0. Conversely if § V(x) dx = and V(x) £ 
0, then this is a sufficient condition for the presence of 
a bound state. This proves the equivalent of Simon's 
theorem for the 4th-order bSDE and, furthermore proves 
the QIC in (3+1) dimensional Minkowski space. 

Note that the version of the QIC that we have proved 
is this: The QFs imply that either ( T^ n (t, 0) ) = 
everywhere along the world-line, or 

<f(T^ n (t,0))^dt>0, (69) 

which we emphasize is slightly stronger than the AWEC. 
Splitting the energy density into its positive and nega- 
tive parts, this implies that as long as ( Tqq"(£, 0) ) is not 
identically zero along the world line, then 

j> < TST(t, 0)) + dx><f( T^ n (t, 0) )_ dx. (70) 

That is, any negative energy "loan" is overcompensated 
elsewhere along the world line. 

V. DISCUSSION 

With the variant of Simon's theorem that we have now 
proved for the bSDE, it is possible to reformulate the 
QIC for a more general range of energy pulses in (3+1) 
dimensional flat space-time, as has been already done for 
the (1+1) dimensional case. In flat space-time, an en- 
ergy pulse which satisfies the QIs (and so the version of 
the QIC discussed above) , must also fulfill an AWEC-like 
inequality. That is, the expectation value of the renor- 
mailzed stress-energy tensor shall violate the condition 
provided by equation ((THJ), and so satisfy (|6^|) . 

Moreover, using the argument derived from equation 
(I2TI1) . the QIC can be extended to (3+1) dimensional 



Minkowski space. Furthermore, from (|2 1 1) and (|70p . we 
know that the positive contributions of the expectation 
value of the renormalized stress-energy tensor must over- 
compesate its negative parts. The key trade-off in the 
current argument is that while we have been able to 
deduce a general result for arbitrary even dimensional 
Minkowski space, and while we are not limited to delta- 
function pulses or pulses of compact support, we have 
on the other hand lost some of the precision informa- 
tion that can be deduced when stronger assumptions are 
made regarding the temporal distribution of the stress 
energy. 

Furthermore, this proof of the QIC that we have given 
in flat four dimensional space-time gives only a partial 
picture of the nature of negative energies and the con- 
straints that can be placed on them, as it does not (yet) 
include the effects of curved spacetime. Indeed, for many 
technical reasons it would be preferable to work with the 
null energy condition (NEC), rather than the weak en- 
ergy condition 18]. That is: A truly complete formu- 
lation of the QIs and QIC should really include curved 
spacetimes, at the very least (3+1) dimensional curved 
spacetimes, and work with some version of the null energy 
condition — this is the key arena wherein the possibil- 
ity of exotic phenomena such as traversable wormholes, 
warp-drives, and time machines are related to the exis- 
tence of negative energies and the constraints thereupon. 
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APPENDIX A: HIGHER DIMENSIONS 

Having now seen the argument in (1+1) and (3+1) 
dimensions it is clear how to generalize to any even num- 
ber of dimensions. Consider a 2m-dimensional (that is, 
([2m— 1] + 1) dimensional) spacetime. Then the QIs are 
equivalent to the statement: 

3 bound state : H = (-1)™-^- + V(x). (Al) 

da;^ m 

Picking the generic class of test functions discussed pre- 
viously, and applying a variational argument, the ground 
state energy is bounded by 

E < j> {V m) (z)] 2 + V(x) [cp(x)} 2 } dx. (A2) 

The kinetic term is 

j[^ m \x)fdx=^ (A3) 
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(where the precise numerical value of k is not important 
as long as it is finite), whereas (assuming all appropriate 
moments exist) 



V(x)\x - n\ n dx. (A4) 



The interchange of the summation and integration is 
again justified by the assumed piecewise analytic nature 
of the test function. Therefore 



2m- 1 



Ea 2m < n 2 + a 2m - 1 ^T x f V(x)\x - n\ n dx 

n=0 ° J 



(A5) 



where there will be some set of constraints on the coeffi- 
cients a n to keep n finite. 

Now since (p(x) is a test function we are always free 
to choose ao 7^ 0, and since probability densities are al- 
ways positive, this forces ao > 0. We are also free to 
(temporarily) choose all the odd £t2 n +i = 0, which is 
a convenience (since it implies the absence of squared 
delta functions coming from differentiations of the abso- 
lute value function) to guarantee k finite. Then 

m— 1 p 

Ea 2m < n 2 + a 2m ~ x ^ x j>V(x)\x - /x| 2 "dz 



n=0 



-0(l/a), 



(A6) 



But then if § V(x) dx < 0, it follows that for a suffi- 
ciently large we can guarantee E <r 2m < 0, whence E < 0, 
and the potential will bind, (thus violating the QIs). This 
is the easy bit. 

• In particular, ^ bound state =>■ <f V{x) dx > 0. 

Now consider the borderline case § V(x) dx = 0. We are 
again free to choose ao ^ (and so ao > 0), and are also 
free to choose both the sign and magnitude of all the am- 
But this implies the potential will bind unless 



In particular 
V/x : 



V(x) \x- fi\ 2n dx = 0; Vrce (0,m- 1). 

(A7) 



V(x) \x-fi\ 



2m-2 



dx = 0. 



(A8) 



Now "turn on" one of the odd a-m+i- Specifically, con- 
sider a% m -i. Since for m > 2 we have 



dx r 



\x\ 2m - 1 ) = Q, 



(A9) 



we can switch on this a2 TO -i coefficient without risk of 
developing a squared delta function in the evaluation of 
k, and so keep k finite. But then since a2m-i is arbitrary 
as to sign and magnitude, to prevent binding we must 
have 



V/i : j> V{x) 



x-{i\ 



dx = 0. 



(A10) 



Differentiating with respect to /i 



V/i : 



V(x) \x-n\ 2m - 2 sign(x-Ai) dx = 0, (All) 



whence, combining the two preceeding equations, we see 

/>oo 

V^: / V(x) \x - y\ 2m - 2 dx = 0. (A12) 
Jn 

Repeated differentiations with respect to fi will now even- 
tually yield V(x) = 0. 

That is, if § V(x) dx = then to prevent a bound state 
we must have V(x) = 0. Conversely, if § V(x) dx = and 
V(x) ^ 0, then there will be a bound state. 

• In particular, ^ bound state 
either V(x) = or § V(x) dx > 0. 



This now is a version of Simon's theorem appropriate 
to the multi-harmonic SDE, leading to a general even- 
dimensional version of the QIC: Either (Tgo"(i,0) ) = 
everywhere along the world-line, or 



:n™( t ,o)), dt>o, 



(A13) 



this last inequality being slightly stronger than the 
AWEC. 
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